Renormalization group analysis of electrons near a Van Hove singularity 
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A model of interacting two dimensional electrons near a Van Hove singularity is studied, using renormalization group 
techniques. In hole doped systems, the chemical potential is found to be pinned near the singularity, when the electron-electron 
interactions are repulsive. The RG treatment of the leading divergences appearing in perturbation theory give rise to marginal 
behavior and anisotropic superconductivity. 



75.10.Jm, 75.10.Lp, 75.30.Ds. 



The relevance of a Van Hove singularity for the high- 
T c compounds was pointed out shortly after their discov- 
ery and a more refined "Van Hove scenario " was 
proposed afterwards [||. In addition there is mounting 
numerical evidence that such a singularity arises in stud- 
ies of strongly correlated 2D systems 0-^. This fact is 
in agreement with photoemission studies in hole doped 
high-T c materials @. 

The dynamics of interacting two dimensional electrons 
near a van Hove singularity are anomalous, as a stan- 
dard perturbative treatment is plagued by logarithmic 
divergences. The study of these divergences requires a 
renormalization group treatment. The feasibility of such 
an approach in condensed matter systems has been an- 
alyzed in ref. jn]]. In the following, we apply a suitably 
generalized version of the scheme proposed in Q to the 
problem at hand. 





FIG. 1. Schematic picture of the Brillouin Zone and the 
Fermi surface of the system considered in the text 

A sketch of the band structure to be discussed is shown 



in fig. 0. 

In a square lattice there are two inequivalent points 
in the Brillouin Zone where a van Hove singularity is 
likely to occur. Near these singularities, the dispersion 
relation is quadratic, and the Fermi surface is well ap- 
proximated by two intersecting straight lines. We will 
not consider the case of perfect nesting, when the Fermi 
surface is straight throughout the entire Brillouin Zone. 
If we restrict our attention to the electronic states near 
these two singularities, we can write a long wavelength, 
effective hamiltonian as: 
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}B (see fig- 



ure |l|). In the following, we will use the parameters 
m = (m X) A + m y) A)/'2 and Am = (m Xy A — m yi A)/2. We 
have omitted spin indices for simplicity. The number of 
possible (isotropic) interactions is restricted by symme- 
try. An umklapp interaction of the type ua^aub^bi is 
allowed. A term like this transfers a pair of electrons 
from the vicinity of one of the singularities to the vicin- 
ity of the other. The total momentum transfer coincides 
with a lattice reciprocal vector. 

The couplings in (|l|)are expressed in units of [ energy x 
area ]. For instance, for a local Hubbard-like interaction, 
U, the parameters in (|l|) are proportional to Ux the area 
of the unit cell. By scaling units in (|l|), all couplings can 
be made dimensionless, that is, u = mu/h 2 . When (|l|) 
is a low energy approximation to a Hubbard like hamil- 
tonian, these dimensionless parameters are proportional 
to -^f, where U is the Hubbard interaction and W is the 
bandwidth. 

We first consider the electrons near a given van Hove 
singularity. An anisotropic scaling of lengths leads to 
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Am = 0. The fact that the interaction is dimensionless 
implies that perturbative corrections should depend log- 
arithmically on the remaining scales implicit in (|l]): the 
high energy cutoff, and the low energy scale, that is, the 
temperature or the external frequency in the diagram. 
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FIG. 2. Diagrams calculated in the text. Full and broken 
lines correspond to propagators of electrons near the two in- 
equivalent singularities. Thin and thick wavy lines denote 
intra- and inter-singularity scattering. 

For instance, the particle hole polarizability is shown 
in fig. (§b). Using a cutoff in energies, such that the 
filled states in the Fermi sea have energies — e c < e < 0, 
we obtain: 
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where e a 



^2/ 2_ 2\ 

Cq — y (note that q is measured from 

the singularity). The angular dependence on q is deter- 
mined by the SO (1,1) symmetry of the dispersion rela- 
tion. Other choices of cutoff, which break this symmetry, 
do not change the low energy behavior of the polarizabil- 
ity, given by ^) . This expression agrees with previous es- 
timates of the polarizability of electrons near a van Hove 
singularity |l^-|T4|]. The imaginary part of the suscep- 
tibility is consistent with the "marginal Fermi liquid " 
hypothesis p5| . 

The inter-singularity polarizability (also shown in fig. 
||b) diverges as: 
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The electron-hole polarizabilities screen the bare in- 
teractions. This flow of the effective couplings makes 
the problem similar to the ID Luttinger liquid [T^] . An 
unrelated renormalization scheme was proposed in [ [l7[ , 
on the basis of the marginal Fermi liquid phenomenology 
[ |l5| . The q dependence of the diagrams derived from ([!]), 
and the electron-electron scattering channels are, how- 
ever, totally different from those proposed in Jl7| ]. 

The intra-singularity effective interaction, to lowest or- 
der, is: 
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This result implies that the system is unstable against 
any distortion with a (small) wavevector along the Fermi 
line. Repulsive interactions favor a spin density wave, 
while attractive interactions lead to a charge density 
wave. Using the same argument for the effective inter- 
singularity interactions, screened by x' •> we a l so find in- 
stabilities for the wavevector (7r, it) which connects the 
two special points. 

To lowest order, we need also consider the renormal- 
ization of the self-energy, given by the Hartree diagram 
shown in fig. (§a). In a conventional Fermi liquid, the 
integrating out of the states in a given energy shell leads 
to a constant contribution, which uniformly shifts the 
one electron levels. The present case differs, in that each 
new energy shell contains more and more states, and the 
shift increases as the renormalization proceeds towards 
the van Hove singularity. For instance, let us define the 
Fermi energy, ep, referred to the singularity. In dimen- 
sionless units, the flow of ep = cf/£c is given by: 
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where I = — log(e c 

The first term in(|5|) is determined by dimensional ar- 
guments. The Fermi energy seems larger when measured 
in terms of the reduced cutoff. The second term gives 
the shift due to the interaction with the high energy part 
of the Fermi sea. We assume that the high energy states 
which are removed lie at energy e c below €f- Unlike in 
normal metals, the second term does not scale in step 
with the first |ll| ]. Thus, for repulsive interactions and 
6p > 0, there may be a non trivial, stable, fixed point, 
at which the flow in ([|) is arrested. As e c is further 
reduced, ep = e* F remains unchanged. In dimensionful 
units, it means that the Fermi energy gets pinned at the 
van Hove singularity. 
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FIG. 3. Schematic view of the Fermi level pinning process 
discussed in the text, a) Density of states in the dopant 
planes, b) Density of states in the CuO planes in the ab- 
sence of the electron-electron repulsion, c) As in b), with the 
interaction treated within mean field, d) As in b), with the 
interaction treated within the RG scheme presented in the 
text. 

A sketch of the physics involved is shown in fig. (||). 
The chemical potential at the CuO planes is fixed by the 
other planes (fig. || a). We know that, for hole doped 
systems, a mean field treatment of the interaction places 
the singularity above the Fermi level (fig. || c). It implies 
that, in the absence of the (repulsive) interaction, the 
singularity should be at lower energies ( fig. || b). Let us 
assume that the singularity is below the Fermi level when 
the interaction is switched off. Then, if the Hartree dia- 
gram is treated using renormalization group techniques, 
starting from the previous situation, the Fermi energy 
will get pinned at the singularity (fig. || d). 

The difference between a mean field and a RG treat- 
ment of the Hartree potential lies in the fact that, in the 
latter case, electrons near the Fermi energy experience 
a larger shift from "fast " electrons, which lie at high 
energies, than from "slow " electrons. The RG analy- 
sis allows for the correlation hole generated by the low 
energy electrons. The pinning discussed here is only pos- 
sible in nominally hole doped systems, in agreement with 
photoemission experiments [ [1Q[ . 

In the absence of a BCS instability (to be discussed 
below), the low order processes that we have considered 
so far renormalize down the bare interactions, in the way 
shown in (||). The hamiltonian ([!]) flows towards a free 
fixed point. In the renormalization process, however, 
other quantities, such as the quasiparticle weight, are also 
modified (wavefunction renormalization). Because of the 
logarithmic divergences in the polarizabilities, the quasi- 
particles are renormalized throughout the Fermi surface. 
The general form of this renormalization is: 
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Thus, the spectral weight of the quasiparticle states 
is energy dependent, and decreases as the Fermi energy 
is approached. This behavior is reminiscent of a ID 
Luttinger liquid 0, and qualitatively agrees with the 
"marginal Ferm liquid " hypothesis p5[ . In addition, the 
transport properties, in the normal state, of the electrons 
near a singularity are strongly modified by the existence 
of a Umklapp electron electron scattering channel jp|,[l8| . 

Finally, we analyze the BCS channels, shown in figure 
(|c) and (gd). These dia grams also induce logarithmic 
corrections to the couplings, which diverge like the loga- 
rithm of the cutoff squared. Thus, this effect dominates 
over the influence of the electron-hole renormalization 
from the diagrams in fig. (||b). 

Because of this anomalous divergence, a straightfor- 
ward approach of RG techniques is not possible. We have 
chosen to scale the couplings as function of I 2 = log 2 (e c ) 
which corresponds to summing the leading log 2 diver- 
gences [||. Defining v as the intra-singularity scattering, 
and v' as the inter-singularity scattering, the flow is given 
by: 
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These equations can be integrated, and the flow follows 
the lines v' = ± (k ± \Jk 2 + v 2 ), where k is a constant. 
The flow of the couplings is schematically shown in fig. 




FIG. 4. Schematic view of the flow of the intra- (v) and 
inter- singularity scattering (v') in the BCS channel. 
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A purely repulsive v is renormalized towards v = 0, 
which is a marginal fixed point. If \v'\ > \v\, the cou- 
plings eventually become large and attractive. The flow 
towards strong coupling is the renormalization group ver- 
sion of the BCS instability. Thus, this instability always 
takes place when \v'\ > \v\, irrespective of the sign of the 
couplings. 

The instability arises because a large v' favors a coher- 
ent superposition of Cooper pairs at both singularities. 
A repulsive interaction ( v' > ) implies that the conden- 
sate amplitude has opposite signs at the two singularities, 
i. e. d-wave pairing. 




a) b) 



FIG. 5. Kohn-Luttinger diagrams which renormalize the 
couplings in the BCS channel. 

Usually, a repulsive interaction in q space is, either 
constant, or decays for large g's. In the present case, 
however, v and v' are screened in very different ways 
by the electron-hole polarizabilities calculated in (H) and 
(||). The relevant diagrams are shown in fig. (||). Intra- 
singularity screening should dominate (provided that 
log (tt^) ~ 0(1)), as can be appreciated by compar- 
ing (J2J) and (||). These processes are the equivalent to 
those originally discussed by Kohn and Luttinger [ [l9| , 
when studying the possibility of superconductivity from 
screened repulsive interactions. The physics is very much 
the same, except that, in the present case, the anisotropy 
of the screening matches well the phase space available 
for the formation of Cooper pairs. Thus, the instability 
arises in a more natural way. In addition, the existence 
of a van Hove singularity modifies the relation between 
the BCS gap and the dimensionless coupling, which now 
becomes A BCS ~ T c ~ e c e~ x l ,{y,v' > 0). 

In summary, we have shown that a system of interact- 



ing electrons near a 2D van Hove singularity can be an- 
alyzed using renormalization group methods. The treat- 
ment presented here is still far from complete, but the de- 
rived features are independent of the details of the model, 
and consistent with the observed behavior of high-T c ma- 
terials. 
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